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Abstract. Inspired by the algorithm for solving the discretisation in time of 
the evolution problem for an implicit standard material, presented in [1], we 
propose a variational formulation in terms of bipotentials. 

1. Introduction 

In the paper [1] Berga and de Saxce propose a bipotential for the constitutive law 
of a soil and further they proceed with the variational formulation of this model. 

We are interested in the precise formulation of the model, especially we want 
to understand from a mathematical viewpoint the recipe proposed in [1] for using 
bipotentials in order to get a variational formulation of their model. We regard 
this as the first step towards the establishment of a general variational theory of 
bipotentials. 

The following paragraph, extracted from [1] page 414, is revealing for two reasons: 
(a) the understanding of the motivation for introducing the bifunctional in order 
to adapt the Uzawa algorithm for implicit constitutive laws; (b) the imprecision 
concerning the understanding of the proposed new algorithm, related to the fact 
that, as we shall see, the simultaneous minimization of the bifunctional is not in 
fact how the algorithm works. 

" One of the advantages of the new formulation is to extend the classical Cal- 
culus of Variations to non associated constitutive laws. In the theoretical frame 
of the Implicit Standard Materials, a new functional, called bifunctional, is intro- 
duced, depending on both the displacement and stress field. The exact solution 
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of the Boundary Value Problem corresponds to the simultaneous minimization of 
the bifunctional, firstly with respect to kinematically admissible displacement fields, 
when the stress field is equal with the exact one, and secondly with respect to stati- 
cally admissible stress fields, when the displacement field is the exact one. The two 
minimization problems are the direct extension of the dual variational principles of 
displacements and stressesJ^ 

The notion of bipotcntial (definition 3.f) has been introduced in [20], in order 
to formulate a large family of non associated constitutive laws in terms of convex 
analysis. The basic idea is explained further in few words. In Mechanics the 
associate constitutive laws are simply relations y £ d4>{x), with 4> : X ^ KU {+00} 
a convex and lower scmicontinuous function. By Fenchel inequality such a relation 
is equivalent with (t>{x) + (t>*{y) = {x,y), where (f>* is the Fenchel conjugate of 4>. 
It has been noticed that often in the mathematical study of problems related to 
associated constitutive laws enters not the function (f>, but the expression 

which we call "separable bipotential" . The idea is then to use as a basic notion 
the one of bipotential b : X x Y ^ M.\J {+00}, which is convex and Isc in each 
argument and satisfies a generalization of the Fenchel inequality. To non associated 
constitutive laws thus corresponds bipotentials which arc not separable. 

There arc many such laws which can be studied with the help of bipotentials, 
as witnessed by the papers listed further. In many of these papers bipotentials are 
used for numerical purposes and several ad hoc algorithms have been suggested 
and exploited for applications. Here is a partial list of constitutive laws which have 
been described by bipotentials: non-associated Driicker-Prager [18] and Cam-Clay 
models [17] in soil mechanics, cyclic Plasticity ([16], [3]) and Viscoplasticity [11] of 
metals with non linear kinematical hardening rule, Lemaitre's damage law [2], the 
coaxial laws ([19], [22]), the Coulomb's friction law [20], [16], [4], [9], [10], [12], [18], 
[21], [13], [7]. A complete survey can be found in [19]. 

Later we started in [5] [6] [7] a mathematical study of bipotentials and their re- 
lation with convex analysis. This paper is another contribution along this subject, 
concerning mathematically sound variational formulations and algorithms for nu- 
merically solving the quasistatic evolution problem for constitutive laws of implicit 
standard materials. For another paper which contains a variational formulation via 
bipotentials for the particular case of separated bipotentials, see [14]. 



2. Notations and prerequisites from convex analysis 

X and Y are topological, locally convex, real vector spaces of dual variables 
X G X and y G Y, with the duality product (■, ■) : X xF — ?> M. We shall suppose that 
X, Y have topologies compatible with the duality product, that is: any continuous 
linear functional on X (resp. Y) has the form x 1— ;> {x,y), for some y €Y (resp. 
y I— {x,y), for some x G X). We use the notations: 

- R = MU {+00}; 

- the domain of a function : X — > R is dam (j) — {x G X : 4>{x) £ M}; 

- To{X) ^{(f):X^R: is Isc and domcj) ^ 0}; 
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for any convex and closed set A C X, its indicator function, ^a, is defined 

by 

if .T e A 

+00 otherwise 



- the subgradicnt of a function : X — > R at a point x ^ X is the (possibly 
empty) set: 

= {u eY I'iz e X {z - x,u) < - (l3{x)} . 

- the inf-convolution of two functions (j),'ip E Tq{X) is the function ^Dt/j G 
To{X) defined by: for any x E X 

(pn\ip(x) = inf {0(it) + t\}{v) : u + v ^ x} 
3. BIPOTENTIALS AND SYNCS 

Definition 3.1. A bipotential is a function b : X x Y ^ M., with the properties: 

(a) for any x € X, ii domb{x, ■) ^ % then 6(x, •) g Tq{X)\ for any y E Y , li 
domb{-,y) ^ then 5(-,y) e ro(y); 

(b) for any x E X,y eY we have &(a;, y) > {x, y); 

(c) for any {x,y) E X xY we have the equivalences: 

(1) yEdb{-,y){x) xEdb{x,-){y) b{x,y) ^ {x,y) . 
The graph of b is 

(2) M{b) = {{x,y)EXxY \ b[x,y) = {x,y)} . 



Bipotentials are related to syncronised convex functions, defined further. 

Definition 3.2. A sync (syncronised convex function) is a function c : X x 
y[0, +00] with the properties: 

(a) for any x E X, \i domc{x, ■) ^ % then c(x, •) E Tq{X)\ for any y E Y , li 
domc{-,y) ^ then c(-,y) E Ta{Y)- 

(b) for any x E X,\i domc{x,-) ^ and the minimum min{c(a;, y) : y E Y} 
exists then this minimum equals 0; for any y E X, ii domc{-,y) ^ and 
the minimum min{c(a;,y) : x E X} exists then this minimum equals 0. 

Proposition 1. A function b : X x Y ^ W is a bipotential if and only if the 
function c : X x Y ^ R, c{x, y) ~ b{x, y) — {x, y) is a sync. 

Remark 1. The string of equivalences (1) justifies the name "syncronised convex 
function", as it expresses the fact that critical points of functions c{x, •) are related 
with critical points of functions c{-,y). 

With the notations from proposition 1, we have AI{b) = c^^(O). Also, for any 
X E X and y eY , property (a) definition 3.2 of syncs is equivalent with: 

epi{c) n {x} xYxR and epi{c) nX x {y} xR 

are closed convex sets, where epi{c) is the epigraph of c: 

epi{c) = {{x,y,r) E X xY xR : c{x,y) < r} 

An interesting fact is that duality products do not enter in the definition of 
syncs. As an application, let {X,Y) be a pair of spaces, (•,•) and (•,•)' be two 
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duality products, defined on X xY, and c: X xY ^ [0, +00] be a sync. We define 
the applications: 

b,b' : X xY ^RU {+00} b{x, y) = c(.t, y) + (x, y) , b'{x, y) = c{x, y) + (x, y)' 

Then 6 is a bipotcntial with respect to (•, •) and b' is a bipotential with respect to 
(•,•)'. As a corollary, if we have a bipotential 6 with respect to the duality product 
(•, •) and (•,•)' is another duality product, then the application b' defined by 

b'{x,y) ^b{x,y) - {x,y) + {x,y)' 

is a bipotential with respect to the duality product (•,•)' and M{b) — M{b') (they 
describe the same law). More generally, we have the following proposition concern- 
ing transformations of syncs. 

Proposition 2. Let (■,■)), (X', F' , (•, •)') be two pairs of spaces with their 

respective duality products, T : X ^ X' and L : Y ^ Y' be two linear, bijective, 
continuous transformations, a > and c' : X' xY' ^ [0, +00] be a sync. Then the 
function 

c : X X y ^ [0, +00] , c{x,y) = ac'{Tx,Ly) 
is a sync and c^^(O) = c'~^(0). 

Proof. The application c' is a sync, therefore it satisfies conditions (a), (b) from 
definition 3.2. It is straightforward that c is convex and Isc in each argument, 
therefore condition (a) definition 3.2 is a consequence of the same condition for c'. 
Also, because T and L are bijective, condition (b) for the application c follows from 
the same condition for c'. □ 

The following is definition 3.1 [5]. 

Definition 3.3. A non empty set M C X x Y is a BB-graph (bi-convex and 
bi-closed) if for any x ^ X and y <EY the sections 

M{x) = {yeY\{x,y)eM} 

M*{y) - {xeX\ {x,y)eM} 

are convex and closed. 

For any BB-graph AI the indicator function '^'m is obviously a sync. To this 
sync corresponds the bipotential 

boo{x,y) = {x,y) + 'i/M{x,y). 

In particular, this shows that to a BB-graph we may associate more than one 
bipotential. Indeed, if M is maximal cyclically monotone then it is the graph of a 
separable bipotential, but also the graph of the bipotential associated to the sync 
(that is a bipotential of the form 600 )■ Therefore maximal cyclically monotone 
graphs admit at least two distinct bipotentials. 

4. Implicit standard materials described by bipotentials 

In the mechanics of standard materials, the evolution problem is generally given 
by a set of equations, inequations, boundary and initial conditions. They can be 
structured in three groups: kinematical equations, equilibrium equations and the 
constitutive law modeling the material behavior. 
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4.1. Notations. The configuration of tlie body is represented by Q,, an open, 
bounded set witli piecewise smootii boundary d^l. 

We denote by n tlic dimension of the configuration space {n = 1,2 or 3), thus 

n c R". 

The boundary decomposes in two disjoint parts: on OqCI displacements are im- 
posed, wliilc given surface forces act on the remaining part of the boundary denoted 
by diil. The closure of f2 is denoted by ft. 

The following quantities are considered. 

- It is the displacement of the body with respect to the configuration J7, 

- e = D{u) = — (Vu + Vm^) is the associated strain. The trace of the strain 

is denoted by e„i ~ — tr e and the strain deviator is 
n 

- The strain e decomposes additively into elastic and plastic strains 



The traces of elastic strain e*^ and plastic strain e'' are denoted respectively 
by e^, ej^, and their deviatoric parts are e^, respectively, 

- The stress field is denoted by tr, its trace is the hydrostatic pressure Sm = 
tr a and s denotes the stress deviator. 

- S is the elasticity tensor modulus. 

- The density of volumic forces is /„ ; on diVt act the surface forces with den- 
sity fs- The class of stress fields a which satisfy the equilibrium equations: 

div (T + /„ = in , a ■ n = f s otl diVL 

is denoted by SA{fy, fs). 

- The imposed boundary displacements on dofl arc denoted by w. In fact, 
it is useful for further computations to consider the imposed boundary 
displacement u to be defined over all f2. The class of displacements u, such 
that u — u = on dfl (possibly in the sense of trace) is denoted by CA{u) 
and called the class of displacements which are kinematically admissible 
with respect to u. 

Let Sym{n) be the space of n x n real symmetric matrices and Sy'mo{n) C 
Sym(n) the subspace of real symmetric matrices with null trace. The decomposition 
of a real symmetric matrix into hydrostatic and deviatoric parts can be expressed 
by the linear bijectivc transformations: 

Tx : Symln) ^ ^ x Symoin) , Ti(e) = ( - tr e,e - -(tr e)Ir^ 

\n n 

T2 : Sym{n) — > M x Symo(n) , TIjfcr) = itra,a (tra)In 

V n 

With the notations previously made, for any strain value £ G Sym{n), or for any 
stress value a € Sym{n), the decompositions in hydrostatic and deviatoric parts 
are: 

Ti(e) = (e„,e) , T2(cr) = (s„i, s) 

(In order to keep track of physical dimensions, we should introduce two spaces 
Sym{n), one for strains and the other for stresses, or introduce units of measure. 
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but we feel that such notations are only making the presentation unnecessary com- 
plicated.) 

We shall consider the following duality products: 

(•, •) : Sym{n) x Sym{n) M , (e, cr) = tr (ecr) 

(•,•)': (M X 5'2/mo(n))x(M X 5'ymo(ri)) , ((e™, e), (s™, s))' ^ e™s™ + (e, s) 

Remark that the first duality product is the one entering in the formulation of the 
dissipation (as an integral over the body configuration D, of {eP,(7. The second 
duality product will be used for the plastic bipotential, see later for the example of 
the Berga & de Saxce bipotential for the non-associative Driicker-Pragcr law. The 
relation between these dualities is: 

(£,a)^(Ti(e),T2(<7))' 

therefore (by passing to associated syncs and back) we can easily transform bipo- 
tentials expressed in coordinates (e, a) into bipotentials expressed in coordinates 

((^mi (-^m; ^))- 

The kinematical equations are: 

(3) e = i (Vu + Vu^) , u e CA{u) 
The equilibrium equations are: 

(4) aeSAif^J,) 

The constitutive equations (besides the additive decomposition of the strain into 
elastic and plastic parts) arc expressed with two bipotentials: the clastic and the 
plastic bipotential respectively. 

The elastic bipotential is defined by the elasticity tensor modulus and it has the 
form: 

(5) &e(e^a) = ^(e^5e^) + i(^-V,a) 

The elastic bipotential is defined over pairs of dual variables (clastic strain, stress). 
It is a separable bipotential, expresses as the sum of the (density of) the elastic 
energy and it's dual. Moreover, this bipotential is quadratic in each variable. 
The plastic bipotential 

(6) bp = bp{iP,a) 

is defined over another pair of dual variables, namely (plastic strain rate, stress). 
In the case of standard materials, the plastic bipotential is separated (expressed as 
the sum of the plastic potential and it's dual). For implicit standard constitutive 
laws which can be expressed by a bipotential (like for example the non-associative 
Driickcr-Pragcr law), the bipotential is not separated. 
The constitutive equations are: 

(7) £ = 

(8) £^ e96e(e^■)(cT) 



(9) 



£^'e96p(£P,-)M 
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The constitutive equation (8) is equivalent with e"^ = S^^a, which is a hnear 
equation. In order to enhance the resemblance between (8) and (9), we could 
differentiate with respect to time in the constitutive equation for e'^ and then express 
the result with the help of the elastic bipotential: 



(10) edbe{i',-){(j) 

5. NON-ASSOCIATED DrUCKER-PrAGER ELASTO-PLASTICITY 

An important example of an implicit standard material is provided by the non- 
associated Driicker-Prager constitutive law. Here we follow the presentation from 
[!]• 

5.1. Plastically admissible stresses. The model is characterized by a Driicker- 
Prager plastic yielding surface. The set of plastically admissible stresses is the 
following cone: 

Kstress = jo" = 7;SmI + s such that -^||.s|| + S,ntg<j) < c 
L 3 kd 

Here c is the cohesion, (j) is the friction angle and kd is a constant whose significance 
is explained in [1] section 3, relations (3.1), (3.2). 

We denote by K'^^^^^^ ~ T{Kstress) the same cone in coordinates (s„i,s) of the 
stresses. 

5.2. Plastically admissible strain rates. Let g [0,0] be the dilatancy angle 
(if 6* = then we are in the case of associated Driicker-Prager elastoplasticity). 
The set of admissible plastic strain rates is the cone: 

Kstrcnn = je^ = \e^rJ + 6^ such that kdtgOWe^ < < 

We denote by A'^train ~ T{Kstrain) the same cone in the representation (6^,76) of 
the strains. 

5.3. The flow rule. The constitutive equation for the evolution of the plastic 
strain has the following expression: 

(11) {{eP, + kd{tg<j> - tg0)\\eP\\),en G (s,„, s) 

Theorems 4.1, 4.2 from [1] are collected into the following. 

Theorem 5.1. Let 6^ : (K x Symo{n)) x (R x Syma{n)) ^ R U {+00} he the 

function: 

(12) 

r CiCn + C2{sm - jf^)\\e\\ if {sm,s) G K'^t^^^^ and 

b'piiem, e), (Sm, S)) = < (e„, e) e K'^trazn 

[ -foo otherwise 
where \\e\\ is the norm defined by ||e|p = (e,e) and the constants 
Ci = -^ , C2 = kditg0-tg(t>) 

are coming from the flow rule (11). Then: 

(a) bp is a bipotential with respect to the duality product (•, •)', 
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(b) the non-associated Driicker-Prager constitutive equation for the evolution 
of the plastic strain (11) can he expressed with the help of the bipotential bp 
as 

&;((eP.,ef),(s„.,s))-((ef„,eP),(s™,s))' 

As an application of proposition 2, we obtain the following characterization of 
the Driicker-Prager constitutive law. 

Corollary 1. In the coordinates (e, cr) G Symin) xSym(n), with the duality product 
{e,a) = tr{ea), the non-associated Driicker-Prager constitutive law (11) can be 
expressed with the help of the bipotential bp : Sym{n) x Sym(n) — >■ M U {+00} 
defined by: 

(13) 6p(e,a) = «'K.,_(a) + «'K,,,„,„(e) + Citr£+C2(ira--^)||e--(ir£)/||- 

tg (p n 

- (^1 - {trs){tra) 

Remark 2. The term containing C2 represents a coupling between the hydrostatic 
part of the stress and the deviatoric part of the strain (rate). If C2 = then 
we get the associated Driicker-Prager constitutive law. In this case the last term 
from the right hand side of the expression (13) can be eliminated by modifying 
the cones K stress and K strain- But if C2 7^ such a modification cannot be made 
because of the coupling between deviatoric and hydrostatic parts, so this last term 
in the expression of bp can not disappear by a modification of the cones K stress and 

^ strain • 

6. Time discretisation of the evolution problem 

Given as the initial data the displacement uq and the initial plastic strain Eq, the 
boundary data u = u(t), fs = fs{t), and the volume forces /„ ~ fv{t), for t G [0, T], 
a solution of the evolution problem is a collection {u,e^,e'^,a) of fields dependent 
on t, which satisfy the kinematical, equilibrium, constitutive equations, as well as 
the initial and boundary conditions. 

We want to give a variational formulation of the time discretisation of the evo- 
lution problem. For this we consider a discretisation 

{io = 0,ti,...,tjv =T} 

of the time interval [0,T]. For each k = 0,...,N we denote by (wfc, e]^, e^, CTfe) the 
unknowns at the moment tk- We shall use also the notation: for any k = 0, ...,iV, 
let Aifc = tk+i — ti;, Aufc = Mfc+i — Mfc, and so on, for all fields, known or unknown. 

Further on, we shall replace the time derivatives from the evolution equation by 
finite differences with respect to the considered time discretisation. The problem 
which we want to solve is the following one. 

6.1. Problem (Pdisc). Given (u^-, e^, e^, Cfc), find (Aw, Ae^, Ae^, Act), solution of 
the following problem: 

(14) Ae" + AeP = D (Am) 



(15) 



Ae"^ = SAa 
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(16) ^A.^ea6,(^A.V)K + A.) 



(17) Aae5^(A/„,fe,A/,,fe) 



(18) Au g CA{Auk) 

The unknowns (uk+i, £^_^_l, £1+1^ '^k+i) arc obtained as 

Uk+i = uk + Au , el^^ = £fe + , ■•■ 

Our first coneern is to express (Pdisc) with the help of bipotentials. 
Lemma 6.1. For any k ~ 0, N — 1, the function 

(19) 6p,fc(A£P, Aa) = Atkbp (^Ae^ , + Aa^ - {Ae^, g^) 

is a bipotential and the equation (16) is equivalent with 

(20) Aef e96p,fe(A£P,-)(Aa) 

Proof. Let us show that 

Cp,fe(A£P, Aa) = fop,fc(A£^ Aa) - (Ae^, Aa) 
is a sync. For this we introduce the sync associated to the bipotential bp, namely 
Cp(Aef, Aa) = 6p(A£f, Aa) - (Ae^, Aa) 

Remark that 

Cp,fe(AeP, Aa) = At^ Cp (^^Ae^, a^ + Aa 

We apply proposition 2 and get the result. By consequence, the function bp_k 
defined by (19) is a bipotential. From here, the second part of the proposition is a 
straightforward computation which is left for the interested reader. □ 

6.2. Simplification of the boundary conditions and volume forces. It is 

not a restriction of generality to suppose that the boundary conditions and volume 
forces are trivial, that is to suppose that equations (17), (18) have the following 
form: 

(21) AaG 5*^(0,0) 



(22) Au G CA{Q) 

Indeed, let us choose a field Aa G SA{Afy^k, Afs^k)- If we define the new 
unknowns: 

Alt' = Au - Au , Aa' = Aa - Aa 

then we could use using again proposition 2 in order to prove that the constitutive 
equations, in the new unknowns, can be expressed by bipotentials. 
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In order not to use a too heavy notation, further on we shaU assume (21), (22) 
and we shah neglect the change of unknowns (thus maintaining the notations Aw, 
Act). 

6.3. Elimination of several unknowns. We can simplify the problem (Pdisc) 
by a standard argument involving the elimination of the unknowns A£'^,Ae'', by 
using an inf-convolution. 

Indeed, let us denote Ae = Ae*^ + Ae^. By equation (14), Ae can be deduced 
from Au. 

For any Act, the functions 6e(-, Act) and 6p_fc(-, Act) are not everywhere infinite, 
are convex and lower scmicontinuous, therefore we can define the inf-convolution 
of them: 

(23) AbkiAe, Act) - {b,{; Act) □ bp^k{; Act)) (Ae) 
Lemma 6.2. 

(24) Act e aAfefe(-, Act)(A£) 

is equivalent with: there are Ae*^, As^ , such that Ae = Ae*^ + Ae^, which satisfy, 
together with Act, the equations (15), 16). 

Proof. Indeed, by a well known property of inf-convolutions, equation (24) is equiv- 
alent with: there are Ae"^, Ae^, such that Ae = Ae*^ + Ae^, which satisfy 

(25) Act G 95e(-,ACT)(Ae'=) 

(26) Aa edbp^k{-,Aa){Aen 

But both be and bp^k are bipotentials, therefore (25) is equivalent with (15) and 
(26) is equivalent with (20), which is equivalent with (16) by lemma 6.1. □ 

Remark 3. Because of the particular form of (quadratic function), the inf- 
convolution Afcfc (• , Act) is differentiable, with Lipschitz gradient, as a kind of Moreau- 
Yosida regularization. Therefore the inclusion (24) is equivalent with a standard 
equality, because the set from the right hand side contains only one element. This 
is an well known advantage of this elimination of unknowns in associated plasticity. 

Let us list the properties of the function Abk' 

- it is lower scmicontinuous (even differentiable, with Lipschitz gradient in 
the first argument) 

- Abk is defined via an inf-convolution of a bipotential of type (13) with the 
elastic bipotential be, therefore it satisfies the same growth inequality as 
be namely there is a constant C > such that for any Ae £ Sym{n) and 
Act G Sym{n), if A6fc(Ae, Act) < +oo then 

A6fc(Ae,ACT) < C (||Aef + HAct^ ) 

where || • || is an arbitrary euclidean norm on the space Sym{n), 

- it satisfies a weak form of the Fenchel inequality, 

A6fc(Ae,ACT) > (Ae,ACT) 

- it is convex in the first argument, but not in the second, therefore it is not 
a bipotential, as it is stated in Theorem 6.1 [1]. Remark however that the 
proof of Lemma 6.2 uses the fact that the function 6p ^ is a bipotential. 
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We collect the partial results obtained so far into the following theorem, which 
provides a simplified form of the problem (Pdisc). 

Theorem 6.3. The problem (Pdisc) is equivalent with the following one: find 
(Au,Acr) e CA{0) X 5A(0,0) which satisfy (24). 

7. VARIATIONAL FORMULATION OF THE PROBLEM (PDISC) 

We give further a variational formulation a la Nayroles [15] of the following 
general problem, which contains (Pdisc) as a particular case. 
We consider a first pair of spaces in duality: 

- X ~ L^{^, Sym{n)) is the space of the deformation fields e, 

- Y = L^(f2, Sym{n)) is the space of stress fields a. 

Instead of equalities X,Y ^ LF'{VL^ Sym{n)), we may consider that X and Y are 
topological, locally convex, real vector spaces of dual variables e £ X and cr G F, 
with the duality product (•, •)! : X x y — > M, endowed with two injective continuous 
linear transformations A : X L^{^, Sym{n)) and B : X ^ L^iP-, Sym{n)) such 
that 



In the integral we see the duality product (scalar product) on the space Sym{n) of 
n X n symmetric real matrices. In the right hand side we see the duality product 



The space X, Y may be finite dimensional (for example associated with a dis- 
cretisation in space by finite elements) or infinite dimensional. In the following 
we shall omit to mention the injections A, B or any other similar transformations 
which may appear. As an exception, in the following theorem 7.1, part (I), we need 
the spaces X, Y to be " large enough" in order to be able to prove that a solution 
of the variational formulation is also a solution (almost everywhere) of the original 
problem. 

U is the space of CA{0) displacement fields u € W'^''^{n,W), with and u = 
on i9o^^ in the sense of trace. The linear transformation D : U X, e = D{u) = 



The space Yq C Y oi statically admissible 5^1(0, 0) stresses appear as the space 
of (7 G y with the property that for any u € ?7 we have 



We consider a function b : Sym{n) x X — 5- R with the following properties: 

(a) K C Sym{n) is a closed convex set of the form K = a + Kq, with a G 
Sym{n) and Kq C Sym{n) a closed convex cone, such that S (this is 
the set of plastically admissible stresses, as in the definition of the Driicker- 
Prager plasticity). Let ttk ■ Sym{n) K he the projection on this cone; 

(b) b is lower semicontinuous in both arguments, differentiable with Lipschitz 
gradient and convex in the first argument; moreover we suppose that the 
Lipschitz constant of the gradient of b in the first argument is continuous 
with respect to the second variable; 

(c) b satisfies, for any e,a G Sym{n), the inequality: b{e, a) > (e, a), 




(scalar product) of with itself. 
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(d) there is a constant C > such that for any e G Sym{n) and cr G A' we have 
(27) bie,a) < C (||ef + ||af ) 

Associated to the function b is the "bifunctional" of Berga and de Saxce: 

B{e,a) = / b{s{x),a{x)) dx 



Our main theorem is the fohowing: 

Theorem 7.1. Suppose that the function b takes only finite values, that is for any 
(£,cr) G Sym{n) x Sym{n) we have b{e,a) < +oo. 

(I) Let u £ U and a G Yq. The pair {u,a) satisfies almost everywhere in 

a G db{-,G){D(u)) 
if and only if for any e £ X we have 
(28) B{D{u),(j) < B{e,a) - {s,a)i 

(II) For any <E U , cr'^ £ Yq there is a sequence {u^^cr^)^ in U x Yq, such that 
for any k G N; 

a. (global condition) for all v E U the displacement u'^'^^ G U satisfies 

B{D{u'^+'),a'^)<B{D{v),a'') 

b. (local condition) the stress a^^^ satisfies almost everywhere in 17 the rela- 
tion 

(7^+^ edb{-,a^){D{u^+^)) 

(III) If a sequence {u^,a^)k from (II) has a subsequence (denoted by same sym- 
bols) such that converges weakly in W^'"^ to u and converges weakly in to 
a, then (w, cr) is a solution of the problem (28). 

Proof. (I) Wc follow the convention: we identify an clement of 17 G L^{^, Symin)) 
(which is an equivalence class of functions) with its representant, defined almost 
everywhere in VL by Lcbesgue theorem. 

Let u £ U and a G Yq, such that we have a G db{-,a){D{u)) almost everywhere 
in ri. Let us take e € X. Then, by integration of the constitutive relation (and by 
the definition of lo), we have 

b{e{x),cr{x)) dx ~ I (e(x), cr(a;)) dx > b{D{u){x),cr{x)) dx 
Jn Jn 

which is exactly the relation (28). 

Conversely, let us start from the last integral inequality, supposed to be true for 
any e G X. Further we suppose also that X — L'^{il, Sym{n)). Let us pick an 
arbitrary xq in the intersection of the Lebesgue sets of D(u) and a. For any open 
ball B{xo,r) C centered in xq we define G X such that Sr = D{u) almost 
everywhere outside B. We obviously get that 



1 



B{xo,r) \ \JB{xo,r) J B(xo,r) 



b{£r{x),a{x)) dx — i {£r{x) ~ D(a){x),a{x)) dx\ > 



(29) >— biDiu)ix),a{x)) dx 

I B{xo,r) 1 JB(xo,r) 
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For any e G Sym{n) we can choose for any r > (but sufficiently small) an Sr such 
that 

r^O I B(.To,r) I 7s(.o,r) 

and such that we can pass to the limit with r to 0, to obtain: 

b{e,a{xa)) ~ {e - D{u){xq),(j{xq)) > b{D{u){xo),<j{xo)) 

This is equivalent with the satisfaction of the constitutive relation almost every- 
where. 

(II) Suppose that for fc G N we defined the element (w'^jCr'^) of the sequence. We 
want to prove the existence of (u'^^^, u'^^^) which satisfy the global condition (a), 
the local condition (b) and cr'^"^^ G lb- 

By the convexity, growth and continuity conditions on b, we easily obtain the 
existence of a minimizer of the functional u E U B{D{u), a^). This proves 
the existence of u^^^ which satisfy the global condition (a). The local condition 
(b) is in fact the definition of a'''^^ . Because of the differentiability and continuity 
conditions on b, it easily follows from a'^ £ Y that (t'"'+^ g Y. We have to prove 
that cr'=+^ G Yq. For this, we choose an arbitrary v £ U and we integrate the local 
condition (b). We obtain that 

B{D{v),a'') - B{Diu''+^),a'') > {D{v) - D{u^+^),a^+^) 

The left hand side of this inequality is non negative (by the global condition) and it 
can be made arbitrarily small, for example by choosing v = u'^'^^ + Aw, for a given, 
but arbitrary w £ U and A > smaller and smaller. As a conclusion we obtain 
that for any w € U we have 

which implies that a'^^^ G Yq. 

(III) Suppose that {vf^ converges, in the given sense, to {u, a). The sequence 
of functionals v t-^ B{D{v), a^) converges in the variational sense to the functional 
V I— >■ B(D{v),a), so, up to the choice of a subsequence, the minimizers of these 
respective functionals (namely the u'^'^^) converge to a minimizer of the latter 
functional. Therefore (w, cr) satisfy the condition 

B{D{v),a) > B{D{u),a) 

for any v £ U . 

The limit cr is in Yq by construction. We can also pass to the limit in the integral 
form of the local condition, which is: for any e £ X 

B{e, a^) - (e - D{u''+^),a^+^) > B{D{u''+^), a'') 

and we get the relation (28). □ 

The previous theorem contains at part (II) an algorithm for finding a solution of 
the problem (Pdisc). This algorithm is the rigorous reformulation of an algorithm 
proposed in [1] section 8. 

However, this theorem can be improved (and will be, in further research) in 
several respects. Firstly, in the case of Driicker-Prager plasticity, the function Abk 
takes also infinite values. In this case the algorithm for solving the problem (Pdisc) 
should take the following form. Let K denote the set of plastically admissible 
stresses. Then: 
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0. initialize {u^,(t^) (for example take them equal to (0,0), 

1. repeat: given {v!',(7^) eU xY, 

a. (global condition) find u*'"'"^ such that for all v G U 

B{D{u^+^),cj'') - {D{u''+^),<7^) < B{D{v),(j^) - {D{v),cj'') 

b. (local condition) define the stress cr*^^^ almost everywhere in Vl by the 
relation 

Wc don't know yet how to prove that such a sequence converges to a solution of 
the problem (28), which is the weak form of problem (Pdisc). 

Secondly, by exploiting the particular expression of the functions b which appear 
in real plasticity problems, we may be able to prove that sequences {u^,(j^) have 
convergent subsequences, for example by a boundedness argument. 

Another, potentially very interesting subject, concerns Coulomb friction. This 
law can be expressed by a bipotential, [20] [7]. It should be interesting to explore the 
corresponding variational formulation, where the bifunctional will contain volume 
integrals as well as surface integrals. Related to this see also the paper [13]. 
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